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Abstract The paper considers the electrodisintegration of the deuteron for kinematic 
conditions of the JLab experiment E-94-019. The calculations have been performed 
within the covariant Bethe-Salpeter approach with a separable kernel of nucleon- 
nucleon interactions. The results have been obtained using the relativistic plane wave 
impulse approximation and compared with experimental data and other models. The 
influence of nucleon electromagnetic form factors has been investigated. 
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1 Introduction 

New experimental data for exclusive electrodisintegration of the deuteron at high mo- 
mentum transfer Q] can be a good instrument for testing proposed relativistic models 
of nucleon-nucleon (NN) interactions. The specific arrangement of the experiment when 
the final state interaction (FSI) effects are minimized allows one to compare results 
of the calculations performed within the plane wave impulse approximation (PWIA). 
Therefore, there is a chance to investigate the influence of nucleon momentum distri- 
butions produced by various models describing observables. 
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During last 15 years several relativistic models of NN interactions have been elab- 
orated [2j-[6]. Generally, they are based on fully relativistic expressions for matrix 
elements depending on four- momenta of the nucleons under consideration. However, in 
most cases there are difficulties caused by the necessity to perform calculations with 
the zeroth component of nucleon relative momentum pq. They are solved by using some 
constraints on po obtained from physical assumptions that limits the applicability of 
these models to low-energy calculations. In practice, they are based on using nonrel- 
ativistic nuclear interaction models (realistic [7KS] or separable potentials [51lll)|). The 
attempt to apply the covariant separable kernel |11II12| fails because of nonintegrable 
singularities which appear when calculations in the high energy region are performed. 
This problem was solved in [T3]-[TS] where the fully relativistic covariant model was 
elaborated. Now it is interesting to investigate how this model describes the electro- 
disintegration of the deuteron at high momentum transfer where relativistic effects are 
assumed to play an important role. 

Another interesting problem is the influence of the used proton and neutron elec- 
tromagnetic form factors. The widely used model is the dipole fit [16]. However, it is 
well known and intensively discussed that the relation of proton charge form factor 
Ge p to the proton magnetic form factor Gm p obtained by the Rosenbluth separation 
technique, differs from the one obtained by the recoil polarization method 17, ISO- To 
describe the results of the latter method, it is necessary to use parametrization [17] for 
Ge p - It should be noted that in this case the Galster parametrization for the neutron 
electric form factor Ge u US] is applied [2]. In this paper we compare the calculations 
with the original dipole fit for the proton and neutron form factors with those where 
the modified Ge p and Geh are used. 

The paper is organized as follows. The formalism to describe NN interactions within 
the Bethe-Salpeter approach with a separable interaction kernel is presented in Sect. 2. 
Sect. 3 considers the calculated cross section. The obtained results are discussed in 
Sect.4. 

2 Formalism 

In the paper the deuteron electrodisintegraton is considered within the Bethe-Salpeter 
(BS) approach [20] with a separable kernel of NN interactions. It is based on the 
solution of the BS equation: 

$ JM (k;K) = ^JU S 2 (k;K) J d 4 P V(k,p- K)$ JM (p; K) (1) 

for the bound state of the neutron-proton (np) system with the total angular mo- 
mentum J and its projection M described by the BS amplitude $ JM . Here the total 
K = k p + k n and the relative k = (k p — k n )/2 momenta are used instead of the proton 
k p and neutron k n momenta. In general, the BS amplitude can be decomposed by the 
partial-wave states through the generalized spherical harmonic y and the radial part 
<j> [12] as: 

$ J a f(k;K {a) ) = J2(yaM(Wc)*f3 Mko, |k|), (2) 
a 

where K(n\ = (Af<j, 0) is the total momentum of the NN system in its rest frame (here 
it is the deuteron rest frame called the laboratory system, LS); is the mass of the 
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deuteron; Uc is the charge conjugation matrix; a, f3 denote matrix indices; a is a short 
notation of the partial-wave state 2S+1 L P j with spin S, orbital L and total J angular 
momenta, p means positive- or negative-energy partial- wave state. S2(k;K) is the free 
two-particle Green function: 

S^fcK) = (i A"-7 + k- 7-m) (1) (i K-j-k-j-m) (2) . 

In these calculations it is more convenient to use the BS vertex function r which 
is connected with the BS amplitude by the following relation: 

$ JM (k;K) = S 2 (k;K)r JM (k;K). (3) 



After using the decomposition of type (01 for the vertex function the relation between 
the $ JM and F JM radial parts can be deduced: 

a(fco,|k|) = ^S o6 (fco,|k|;a)g 6 (fc ,|k|), (4) 



where S a b is the one-nucleon propagator |12j . To solve the BS equation (JXJ) , we have 
used the separable ansatz for the interaction kernel 

N 

V o6 ( Po ,|p|;fco,|k|;a)= ]T \ij(s)g [ ° ] ( P o, \p\)gf ] (k , |k|), (5) 

»,J=1 

where JV is a rank of the kernel, are model functions; A is a parameter matrix 
satisfying the symmetry property Ajj-(s) = \ji(s); k \p] is the relative momentum of 
the initial [final] nucleons; s = (p p + p n ) 2 where p p is the outgoing proton and p n is 
the neutron momentum, respectively. If the radial part of the vertex function r ,JM is 
written in the following form: 

N 

ffa(po,|p|)= Yl A y( s )s! Ql (Po,|p|)c 3 (s), (6) 

i>i=l 

the initial integral BS equation JT} is transformed into a system of linear homogeneous 
equations for the coefficients Ci(s): 

N 

<=*(«)- Yl /l ifc(s)A fci (s)c i (s) = 0, (7) 
k,j=l 

where 

fl] a] (* ,|k|) ff t a] (fc ,|k|) 



hij(s) = 



i f „ f , 2 ,„ , 9 i (ko, k 5 J (fc , k ) 

--r^r> dfco k dk 2 — (g 

4?r 3Z^y u / (V5/2 - ^ + ie) 2 - *2 



and i^k = Vk 2 + m 2 . Using (|4| and taking into account only positive-energy partial- 
wave states for the deuteron 3 S^~, , the radial part of the BS amplitude can be 
written as follows: 

rh lb ith - gq( fc o ; |k|) /q\ 

MkM -(M d /2-E k + l er-kl (9) 

Thus, using separable g functions we can calculate observables describing the np sys- 
tem. 
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3 Cross section 

When all particles are unpolarized the exclusive d(e, e'n)p process can be described by 
the cross section in LS: 

d 3 cr CMottTn 



dQ 2 d\p n \dn n ~ 2(2tt) 3 M d E e E' e 

x [looWoo + l++(W++ + W-) + l+- cos 20 2ReW+- 
sin 2(f> 2ImW+- - lo + cos 2Re(W 0+ - W -) 
-4] + sin0 2Im(W + + W -)] , (10) 

where CMott = (a cos |/2-E e sin 2 lj) 2 is the Mott cross section, a = e 2 /4ir is the fine 
structure constant; E e [E e ] is the energy of the initial [final] electron; Q' e is the outgoing 
electron solid angle; 6 is the electron scattering angle; Q 2 = —q 2 = —u 2 + q 2 , where 
q = (uj, q) is the momentum transfer. The outgoing neutron is described by momentum 
Pn and solid angle fi n — (9n,4>) with zenithal angle 6 n between q and p n momenta 
and azimuthal angle 4> between the (ee') and (qp n ) planes. The photon density matrix 
elements have the following form: 



,0 Q 2 ,o Q [O 2 



q |q|v2 V q 2 

,o , 2 6 Q 2 o Q 2 f -.,\ 

*++= tan 2 + 2^' l+ - = ~2tf- 

The hadron density matrix elements 

W xx , = W^e x , u , (12) 

where A, A' are photon helicity components |21| . can be calculated using the photon 
polarization vectors e and Cartesian components of hadron tensor 

WW = | Y, \<np:SM s \j»\d:lM>\ 2 , (13) 



where S is the spin of the np pair and Mg is its projection. The hadron current in 
(|13|) can be written according to the Mandelstam technique [22] and has the following 
form: 

< np : SM s \j„\d : IM >= i £ f ^P S p{A(C^ SMs (p CM -P CM )A(£.) 

xrW( fl) s<-> - (-i)*p- 1) r M (p + (-ir|;^ ( o))} (M) 

within the relativistic impulse approximation. The sum over r = 1,2 corresponds to 
the interaction of the virtual photon with the proton and with the neutron in the 
deuteron, respectively. Total P CM and relative p M momenta of the outgoing nucleons 
are considered in the final np pair rest frame (center-of-mass system, CM) and can be 
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written in LS using the Lorenz-boost transformation along the q direction. The Lorenz 
transformation of np pair wave function ipsM s from CM to LS is: 

where r\ = q 2 /s. The interaction vertex is chosen in the on-mass-shell form: 

/>(<?) = 7 „-Fi(9 2 ) - ^ but ~ flu) F2(q 2 ), (16) 

here F\(q 2 ) is the Dirac form factor, ^(g 2 ) - Pauli form factor. The form factors are 
described by the dipole fit model [IB] or modified dipole fit [171119"] . If the outgoing 
nucleons are supposed to be non-interacting then this is the so-called plane-wave ap- 
proximation. In this case the np pair wave function can be written in the following 
form: 

^sm s (p; p) -> 4m s (p.p*; p ) = &*fxSMs (Pi p)S{p - P*), (17) 

where p* — (0, p*) is the relative momentum of on-mass-shell nucleons, XSM S describes 
spinor states of the pair. Taking into account representation (|17[) . the hadron current 
(|14|l can be transformed into a sum: 

< np : SMsW : IM >= i {MC-^XSMs (p* CM ;P CM ) A(C)r^ r) (q) 

r=l,2 

^ - (-i)V - f) r M (p* + (-Iff; K (0) ) } . (18) 

In this paper the cross section of the exclusive electrodisintegration of the deuteron 
d 2 a I dQ 2 d\p n \ [I] is calculated. It can be obtained from (|10|l after integration over the 
neutron solid angle: 

d 2 a f d 3 a 

dQid\p n \ J dQ2d\ Pn \dn n ai ' n - 

On 

According to [T] the integration is performed over Q n : 20° ^ 6 n ^ 160° , 0° ^ <f> ^ 360° . 
Four different Q 2 are considered. The obtained results are discussed in the next section. 



4 Results and discussion 

In this paper the exclusive cross section of electrodisintegration (|19p for kinematic 
conditions of the JLab experiment [T] has been calculated within the Bethe-Salpeter 
approach with the rank-six separable kernel MY6 [15] . The calculations have been 
performed within the relativistic PWIA. The obtained results have been compared 
with experimental data and two theoretical models, the nonrelativistic Graz II (NR) 
[23] and relativistic Graz II QT] separable interaction kernels. 

Figs. I H4l illustrate the cross section depending on outgoing neutron momentum p n 
for Q 2 =2, 3, 4, 5GeV 2 , respectively. The dipole fit model [16] for the nucleon elec- 
tromagnetic form factors has been used. One nonrelativistic Graz II (NR) and two 
relativistic MY6, Graz II separable kernels of NN interactions have been investigated. 
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A good agreement with the experimental data can be seen at low neutron momenta 
|p n | <0.25GeV/c on the figures. The discrepancy between the theoretical models and 
the experimental data increases with jp n j >0.25GeV/c for all the considered mod- 
els. However, we see the agreement of the relativistic models (MY6, Graz II) with 
the experimental data at high neutron momenta. Moreover, the relativistic description 
becomes better with Q 2 increasing and theoretical curves go practically along experi- 
mental points at Q 2 = 5GeV 2 . Therefore, relativistic effects play an important role in 
description of the deuteron electrodisintegration at high momentum transfer and high 
neutron momenta. 

The calculations with the modified proton Ge p [IZ] an d neutron Ge u [H] form 
factors at Q 2 =2 GeV 2 (Fig. [5} and Q 2 =5 GeV 2 (Fig. [6]) are compared to those obtained 
using the dipole fit model for nucleon form factors. Two relativistic models of NN 
interactions MY6 and Graz II have been considered. All the theoretical calculations 
agree with the experiment at |p„| < 0.25GeV/c and begin to deviate from it with 
p n increasing. We can also see slight difference between the cross sections obtained 
using the dipole fit and modified dipole fit models for the nucleon electromagnetic 
form factors. It is interesting that the results calculated within the dipole fit model, 
which does not describe the behavior of the electric form factor of the proton at high 
Q 2 , are virtually undistinguishable from those obtained with the modified Ge p [17] . 
However, the final conclusion which model is better can be made only when the final 
state interactions, negative-energy partial-wave states (P waves) and two-body currents 
(TBC) are taken into account. 

It should be noted that the behavior of the calculated cross section is similar to the 
behavior of the corresponding wave function for the deuteron 3 S^~ partial- wave state 
which is shown in Fig(7] It is seen from the comparison of the cross sections calculated 
with (MY6, Graz II) and without 3 D^ partial-wave state (MY6-S, Graz II-S) in the 
deuteron (Fig|8} that the influence of the 3 state is maximum at low and medium 
neutron momenta. The minima of the 3 S? wave functions are noticeable in the cross 
section. They are smoothed by the 3 D^~ state when the both partial- wave states are 
taken into account. The role of the state increases at high p n . 

It is seen from Figs[T]|5]that the Graz II and MY6 models give qualitatively the sim- 
ilar description of the experimental data within the used approximation. The difference 
between the model calculations and the experimental data can be probably eliminated 
when FSI, P waves and TBC are taken into account. It should be emphasized that 
there is no possibility to improve a theoretical description using the Graz II model. 
As it was mentioned above, FSI is impossible to calculate with the Graz II kernel for 
high-energy particles unless the po component is constrained by some assumption like, 
for instance, in quasipotential approaches. On the contrary, it is possible to take FSI 
into account using the MY6 kernel without constraining pq. 

We should also comment the results obtained using the Paris potential model [S] 
in [I]. A good agreement with the experimental data was achieved when the FSI and 
meson-exchange current effects were taken into account. However, it seems questionable 
to use the nonrelativistic potential elaborated to describe the np elastic scattering 
data for laboratory energies of the colliding particles less than 350 MeV when the 
cross section at high Q 2 and |p„| is calculated. The MY6 model has two important 
advantages. Firstly, it is fitted to describe all available elastic np scattering data [15] . 
Secondly, it allows one to perform calculations without any necessity to constrain po 
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Fig. 1 The cross section (1 1 9 H depending on neutron momentum p n is considered for Q 2 = 
2±0.25GeV 2 . Calculations with the Graz II (NR) [23] (purple dash-dot-dotted line), Graz II 
1 1 1 1 (brown dash-dotted line) and MY6 1151 (red solid line) models are present. The dipole fit 
model 1161 for nucleon form factors is used. 




p n (GeV/c) 

Fig. 2 As in FigJTJ but for Q 2 = 3 ± 0.5 GeV 2 . 

In this paper the comparison of three different models of NN interactions has 
demonstrated that relativistic effects play an important role in the description of the 
deuteron electrodisintegration at high momentum transfer. The result is slightly de- 
pendent of the model used for the proton and neutron electromagnetic form factors. 
Further investigation is required to conclude which models of NN interactions and 
nucleon electromagnetic form factors are reasonable. In particular, it is necessary to 
calculate FSI, P waves and so on. 
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